Magnetic exchange mechanism for electronic gap opening in graphene 
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We show within a local self-consistent mean-field treatment that a random distribution of magnetic 
adatoms can open a robust gap in the electronic spectrum of graphene. The electronic gap results 
from the localization of the charge carriers that arises from the interplay between the graphene 
sublattice structure and the exchange interaction between the adatoms. The size of the gap depends 
on the strength of the exchange interaction between carriers and localized spins and can be controlled 
by both temperature and external magnetic field. Furthermore, we show that an external magnetic 
field creates an imbalance of spin-up and spin-down carriers at the Fermi level, making doped 
graphene suitable for spin injection and other spintronic applications. 

PACS numbers: 71.27.+a,73.20.Hb,75.30.Hx 



Graphene is a two-dimensional sheet of carbon atoms 
whose elementary excitations are massless Dirac fermions 
Due to the vanishing density of states (DOS) near 
the Dirac points, adatoms, such as H, which arc other- 
wise non-magnetic and bond strongly to carbon atoms, 
can easily form local magnetic moments in graphene Q. 
The combination of properties such as ballistic transport, 
extremely large spin coherence lengths 0] , and gating ef- 
fects makes this material very attractive for spintronic 
applications Q. 

For certain applications, such as in the fabrication of 
transistors, it would be desirable to open a large gap in 
graphene. However, opening a gap in suspended sam- 
ples has turned out to be unexpectedly difficult. Dirac 
fermions have long inelastic mean free paths and are typ- 
ically insensitive to disorder and localization effects. The 
Dirac points seem also very robust against instabilities 
in general. Localization effects, nevertheless, can be also 
triggered by scattering mechanisms that operate in the 
spin channel, which is not very sensitive to the relativis- 
ts nature of the electrons. Due to the presence of sub- 
lattice quantum numbers, when magnetic adatoms are 
adsorbed on top of graphene, the interaction between 
two local spins siting on opposite sublattices is antiferro- 
magnctic, whereas spins on the same sublattice interact 
ferromagnetically [a-Q. The interplay of magnetic in- 
teractions and disorder cangive rise to the formation of 
spin textures and clusters [9|, which can localize charge 
carriers as a result of the exchange interaction between 
local and itinerant spins. In this Letter, we show that 
the adsorption of magnetic adatoms constitutes a sur- 
prisingly flexible way to localize charge carriers and open 
a controllable gap in graphene. 

We address this issue by considering a tight-binding 
model, in which local magnetic moments are randomly 
distributed on a honeycomb lattice, and interact with the 
spin of the itinerant electrons through a local Heisen- 



berg exchange interaction; this effective interaction re- 
sults from the interplay between hybridization effects and 
Coulomb interactions, similarly to the Kondo Hamilto- 
nian for metals [Toj | . We assume a self-consistent semi- 
classical picture where the local spins are treated clas- 
sically but the itinerant electrons are treated quantum 
mechanically, accounting for the realization of the dis- 



order in the localized spins |12(. In the numerical 
results, we show the emergence of a large localization 
gap that is driven by the concentration of adatoms and 
the exchange coupling, and which can be controlled by 
temperature and external magnetic field. Furthermore, 
when the chemical potential is shifted above or below 
the gap, we propose that the application of an external 
magnetic field can spin polarize the charge carriers at the 
Fermi level, making doped graphene useful for spintron- 
ics applications. Although this effective description does 
not account for the change of the energy spectrum due 
to the covalcnt character of the bonding and lattice dis- 
tortions caused by the incorporation of the adatoms [HI , 
this model is justified when the magnetic exchange is the 
dominant interaction. In graphene, the exchange cou- 
pling of the Kondo Hamiltonian has been predicted to 
be as large as a few eV [H[ . 

We start from a random distribution of magnetic im- 
purities (or spins) sitting on top of carbon atoms. For 
spins on top carbon sites, the exchange Hamiltonian in 
graphene has a simple form [ljj], and the total Hamilto- 
nian, including charge carriers and external fields, is 



(1) 



where a,i a (bi a ) annihilates an electron with spin a on 
site i of sublattice a = A (B). The first term in Eq. ([T]) 
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describes charge carriers hopping between nearest neigh- 
bor graphene sites, with t ~ 2.7 eV while the second one 
couples localized (Si) and itinerant (s*j) spins through the 
local exchange constant Ji > ; the last two terms in 
Eq. ([1]) describe the interaction of the carriers and spins 
with an external in-plane magnetic field H. Although 
the exchange interaction is local, the indirect interaction 
among the localized spins is long ranged, which, in prin- 
ciple, allows the formation of a true ordered state in two 
dimensions, in spite of thermal fluctuations. 

At the mean-field level, Si ■ Si -t S(i)s* + S-s(i) - 
S(i)s(i), where the magnetizations S(i) = (Sf) and 
s(i) = (sf) must be computed self-consistently at each 
individual site jll|, [l2j ■ The mean-field decomposition 
results in an effective Hamiltonian of charge carriers that 
can be fully diagonalized for any configuration of local- 
ized spins. We then calculate the expectation value of 
the carrier spins s(j) on all lattice sites. The expec- 
tation value of each local spin can be computed back 
by enforcing S(i) to be equal to the local magnetiza- 
tion due to charge carriers at each impurity site, in the 
presence of the external field, S(i) — Bs(0Hi), where 
Hi = J(sf) — Q^-bH and Bg(x) is the Brillouin function. 
The process is iterated until self-consistency is achieved 
for all localized spins (llj . As a result, for a given random 
configuration of localized spins (with concentration n s ) 
and fixed density of carriers n e , we determine the chem- 
ical potential, the energy, E nai on the ro-th level and 
the corresponding probability amplitude, |f/w,o-0)| 2 i for 
occupation of a carrier state with spin uona given site j 
on sublattice a. With these values, we extract the mag- 
netizations of the localized and itinerant spins on each 
sublattice, M a and m a , respectively, and the DOS. The 
classical spins S(i) range from to an arbitrary value 
S m ax- Here we set S ma x = 1/2. 

We have performed calculations for systems with up to 
N=2 x 65 x 65 sites, and used random distributions of 
magnetic impurities corresponding to a given spin con- 
centration n s ; from now on, all quantities should be un- 
derstood as configurational averages taken in ensembles 
varying between 10 and 50 disorder realizations, which 
suffice to yield small enough error bars. As we are in- 
terested in the case where the Fermi level is close to the 
Dirac point, in the following we set n e = 1 (one electron 
per carbon). 

As a check of our scheme, we verified that calculations 
for full coverage of local moments (n s = 1) at fx = 
correctly reproduce the AFM ground state predicted by 
perturbative methods [1,0]; at low temperatures, our re- 
sults also agree with those obtained through Monte-Carlo 
simulations with frozen charge degrees of freedom Q . For 
a random distribution of impurities, Fig. [Ha) shows one 
localized spin configuration for n s < 1 and = 0; we 
see that both FM and AFM clusters are formed. The 
corresponding carrier DOS is depicted in Figfljb): a gap 
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FIG. 1: (Color online) (a) Typical random spin configuration 
(spin up in black, down in orange) for n s = 0.4, N = 2 x 
65 x 65: (/i = 0); (b) density of states (DOS) averaged over 
disorder realizations represented by panel (a). The inset in 
panel (b) shows the DOS over a wider energy range. Vertical 
lines indicate the Fermi level. All data are for J = t, and 
k B T = 0.002f. 



clearly opens around the Dirac point [l6, 17|. 

The mechanism of gap opening in this material 
strongly resembles the double-exchange mechanism in 
manganites [H| . Hopping does not flip the electron spin 
while it moves from one sublattice to another; since 
neighboring adatoms (occupying different sublattices) 
are antifcrromagnetically aligned, it costs an energy J 
for an electron to move between sublattices, leading to 
electronic localization. For isolated adatoms, the aver- 
age energy cost, and hence the localization length, will 
scale with the dilution. In Fig. Gl^a), we show the in- 
verse participation ratio (IPR), which measures the de- 
gree of localization of states, and is defined by the quan- 
tity P an = Y^f a |Vw<x(«)| 4 - For an extended state, P an 
scales with 1/N, whereas for a true localized state it does 
not scale with the system size, and its IPR is compara- 
tively larger. States near the gap have a significantly 
larger IPR than those with energies away from it. There- 
fore, the former states are localized, while the latter, 
whose IPR clearly scales with N, are extended. In the 
dilute limit, these localized states resemble self-trapped 
magnetic polarons, similar to the ones proposed by Na- 
gaev (lij : carriers are localized by exchange interaction 
in spin-polarized droplets located at energetically favor- 
able positions, whose sizes scale with 1/J (sec Fig[2](b)- 
(c) ) . Our numerics also indicates that the localization of 
the charge carriers is robust against local gating effects 
induced by the local reconstruction of the electronic den- 
sity at the impurity sites. We considered this effect by 
adding a diagonal term in the graphene Hamiltonian that 
changes the on-site energy at sites where the magnetic 
ions are attached. We find that the gap is still present 
for on-site energies < J. Since the exchange coupling 
can be extremely large [13, 15 1, the localization effects 
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FIG. 2: (Color online) (a) Inverse participation ratio (IPR) 
calculated under the same conditions as in Fig 1(a) for three 
system sizes: N = 4050 (green triangles), N = 6050 (red 
squares) and N = 8450 (black circles); (b)-(c) Typical eigen- 
state with energy in the vicinity of the gap for J = t (b) and 
J — 2t (c). The size of the symbol is proportional to the 
relative weight of the eigenstate in a given site (spin up in 
black, down in orange). Dark and orange regions represent 
spin polarized droplets (see text). 



discussed here may be dominant for a variety of adatoms 
in graphene. 

We first examine the behavior of the staggered magne- 
tizations, M s = \Ma — Mb\/S and m s = \tjia — msl/s, 
corresponding to local and itinerant spins, respectively. 
From Fig. EJa), we see that due to the nature of the 
RKKY interaction in graphene, even for a diluted system 
with ferromagnetic clusters [see, e.g., FigQJa)], M s at- 
tains its maximum value, n s , at low temperatures. The 
spins of carriers and magnetic impurities are antifcrro- 
magnetically coupled to each other by the exchange in- 
teraction, and hence both magnetizations vanish at the 
same critical temperature T c , for a given concentration 
of adatoms (n s ). Figure [^b) shows the temperature de- 
pendence of the gap, A(T), from which the correlation 
between the curves in panels (a) and (b) is clear: A de- 
creases from its saturated value (proportional to n s J [l6[ 
) at low T, tracking the behavior of both M s and m s , un- 
til it reaches a minimum value when M s = m s ~ 0; due 
to the finiteness of the system, this minimum value is not 
zero, hence setting a minimum energy scale. Below T c , 
the resistivity of the system displays an activated behav- 
ior, typical of a semiconductor for ksT <C A(T = 0), but 
undergoes a metal-insulator (MI) transition at interme- 
diate temperatures. The strong temperature dependence 
of the gap and the temperature-induced MI transition 
resemble the transport characteristics of Kondo insula- 
tors m. 

When a sufficiently strong magnetic field H c is ap- 
plied, aligning the local spins, electronic hopping is no 
longer hindered: the charge carriers dclocalize, and the 



system becomes metallic. This behavior can be seen from 
Figs[3]Jc) and (d), which illustrate that for a fixed T, the 
gap also tracks the behavior of M s with H: A vanishes 
when the spins in both sublattices point in the same di- 
rection. This dependence has important consequences for 
the transport properties of magnetically doped graphene: 
the system will display a very strong negative magnetore- 
sistance, since H tends to close the gap and to increase 
dramatically the conductivity of the system. This ef- 
fect is very similar to the one observed in manganites 
and in magnetic semiconductors, such as EuSe [1 81 ] and 
magnetic hexaborides [2l[; it could also have bearings 
on the colossal negative magnetoresistance observed in 
dilute fluorinated graphene [221 ] . 

As it can be seen from Fig. [3l magnetic field, tem- 
perature, and concentration of magnetic adatoms can be 
used to control the magnitude of the localization gap in 
graphene. In order to manipulate A, it is necessary to 
establish the range of parameters for which the gap opens 
or closes. To this end, we have varied the adatom con- 
centration, n s , and determined the critical temperature, 
T c , and the critical field, H Cl for fixed temperature. The 
critical lines in the main panel of FiglHJa) scale quadrat- 
ically with J, T c oc J 2 , for n s = 1, as seen from the inset 
of Fig. Hta). As one increases the adatom concentration 
from n s = 0.2, the carrier polarization tends to increase, 
which, in turn, enhances the correlation among the local 
moments. As expected, T c and H c rise monotonically 
with the amount of disorder [see Fig. H{b)]. 

In the numerics, we find that both the energy gap and 
the DOS are not strongly dependent of fi. Since the exter- 
nal magnetic field breaks the degeneracy between spin-up 
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FIG. 3: (Color online) (a) Staggered magnetization versus 
temperature for localized spins and carriers and different im- 
purity concentrations: n 3 = 1 (black squares for M 3 , and full 
line with no symbols for m s ), n a — 0.5 (red circles and dashed 
line) and n 3 = 0.3 (blue triangles and dotted line); (b) Gap 
(A) as a function of the temperature [the color code for n 3 in 
(a) is the same in all panels]; (c) Staggered magnetization as 
a function of applied magnetic field H ; (d) Gap as a function 
of H. J — 0.5t, fi = and N = 2 x 50 x 50 in all cases. 
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and spin-down states, a controlled variation of fx, caused 
e.g., by a gate voltage, can drive the system from an 
insulator to a metallic regime with spin-polarized carri- 
ers, where only spin-up (or spin-down) states are avail- 
able. This is illustrated in Fig.[5](a), where we depict the 
DOS at finite magnetic field for fj, = (dashed vertical 
line), when the system is an insulator, and for \i ~ 0.2i 
(solid vertical line) , when it becomes metallic with a spin 
polarized density of charge carriers at the Fermi level. 
Hence, the application of a external gate voltage, Vg, 
can be used to change the spin polarization of the states 
at the Fermi level, permitting the generation and control 
of spin-polarized currents through the use of randomly 
distributed adatoms in graphenc. 

Noting that a single graphene layer can be separated 
in different regions, with and without the presence of 
adatoms, if during the deposition process part of the layer 
is covered with a mask, one can produce and test spin 
polarized currents using a device sketched in Fig. [5jb) . 
Half of the layer is covered with magnetic adatoms, with 
which a non-magnetic (normal) contact is made, while 
the clean part has a magnetic contact, such as Cobalt. 
An external field is turned on and produces the spin po- 
larized shift in DOS, while the gate voltage moves /j, with 
respect to the Dirac point [Fig. Ufa)]. If a current passes 
through the system from the normal to the magnetic con- 
tact, spin-polarized currents are generated in the magnet- 
ically doped part of graphene and injected in the clean 
part, where it will be tested by the magnetic contact. Dif- 
ferently from other spin-injection schemes experimentally 
tested in this case injection has no external agents, 
thus eliminating unwanted sources of scattering such as 
interface mismatch between graphene and other materi- 
als, like magnetic metals. 

These ideas can also be used to propose a spintronics 
device similar to a spin- valve [2^]. Instead of separat- 
ing the layer in two parts, as above, it is now separated 
in three: two magnetically doped parts separated by a 
clean one. The doped parts are subject to independent 
Vg's and both have normal metallic contacts. With the 
same H and different Vg's, each magnetic part can have 
a specific spin imbalance at the Fermi Level. If a current 
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FIG. 4: (a) Critical temperature T c (M s = 0) versus spin 
concentration n s . (b) Critical field T c versus n a for T — 
O.OOlt. J = 0.5t N = 2 x 36 x 36 The inset in panel (a) 
shows T c in function of J. 
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FIG. 5: (a) Electronic density of states at finite magnetic 
field for J = t,n s = 0.5, [iH = 0.12f, k B T = O.OOli. Dashed 
vertical line: /i = 0, corresponding to an insulating state; 
solid vertical line: /j, = 0.2t (metallic). Black (dark) curve 
represents spin up and orange (light) curve, spin down, (b) 
A spin injection device (see text). 



passes through the system in both directions, the magne- 
toresistance can be controlled by the Vg's . The device 
also produces spin-polarized carriers in one side, and tests 
the polarization in the other. In the device sketched in 
Fig. 5 (b), for a spin concentration of n s = 0.1, and 
exchange coupling as low as J = 0.1 eV, spin polarized 
carriers can be produced by fields less than IT, combined 
with gate voltages of the order of just a few mV. 

In conclusion, we have established that a random dis- 
tribution of magnetic adatoms can localize charge carri- 
ers and open a gap in graphene. The mechanism of gap 
opening is the hindering of hopping due to the antiferro- 
magnetic correlation of magnetic impurities on opposite 
sublattices. The size of the gap depends on the concen- 
tration of adatoms, and on the strength of the exchange 
coupling; it can therefore be controlled by temperature 
and external magnetic field. These two adjustable pa- 
rameters can drive graphene into a metal-insulator phase 
transition in the presence of magnetic disorder. Further- 
more, we showed that a magnetic field can be used to 
produce spin polarized carriers, suitable for spin injec- 
tion and other spintronic applications. 
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